Low temperature thermodynamics of charged bosons in a random potential and the 

specific heat of La2-xSrxCuOi below Tc. 
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We propose a simple analytical form of the partition func- 
tion for charged bosons localised in a random potential and 
derive the consequent thermodynamics below the superfluid 
transition temperature. In the low temperature limit, the 
specific heat C depends on the localisation length exponent 
u: C is linear for v < 1, but for i/ > 1 we obtain C oc T^'". 
This unusual sub-linear temperature dependence of the spe- 
cific heat has recently been observed in La2-xSrxCu04 below 
Tc. 
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I. INTRODUCTION 

A charged 2e Bose liquid of small bipolarons -as a 
microscopic model of the ground state of cupratesEl ex- 
plains the high value of Tc and its 'boomeraiig' doping 
dependencea, the divergent upper critical fieldEI, the nor- 
mal state in-plane and 'c'-axis transport, the magnetic 
susceptibility and the pseudogapQ. The (bi)polaronic na- 
ture of carriers in these materials is supported by obser- 
vations of the characteristic polaronic spectral function in 
the infrared spectraQ, of the isotope effect on the carrier 
masaj, gf flat bands in the high-resolution photoemission 
spectraQ, and by a few other experiments sueh as the 
pair-distribution analysis of neutron scatteringH. 

The bipolaron model overcomes a fundamental prob- 
lem with any theory of high- Tc which is the short-cange 
Coulomb repulsion. As stressed by Emery et ala the 
theories involving real-space pairs with a pure electronic 
(exchange) mechanism of pairing are a priori implau- 
sible due to the strong short-range Coulomb repulsion 
between two carriers. The direct repulsion is usually 
much stronger than any exchange interaction. The at- 
tractive potential generated by the electron-phonon in- 
teraction of the Holstein model may overcome the short- 
range Coulomb repulsion, but inevitably involves a huge 
carrier mass enhancement (otherwise the phonon fre- 
quency would be-jextremely high). On the other hand, 
one of us showedH that the Frohlich electron-phonon in- 
teraction can provide intrinsically mobile intersite small 
bipolarons which are condensed at high Tc (of the order 
of lOOi^). This interaction, operating on a scale of the 
order of leV, compensates the Coulomb repulsion and 
allows the deformation potential (together with an ex- 
change interaction of any origin) to bind two holes into an 
intersite mobile bipolaron in the Cu02 plane. The mass 



renormalisation appears to be smaller by several orders 
of magnitude than in the Holstein model (with the same 
attractive potential). Hence the charge 2e Bose liquid, 
being already an important reference system of quantum 
statistics, is now of particular physical interest. Because 
doped Mott insulators are disordered, the localisation of 
carriers by a random potential plays an important rolej-in 
their low-temperature thermodynamics and transporta. 
In this paper we study the thermodynamics of charged 
bosons in the superfluid phase assuming that some of 
them are localised by disorder. The fact that in the su- 
perconducting phase the chemical potential /x is pinned 
at the mobility edge, /i — Ec, significantly simplifies 
the problem. By using a plausible form of the single 
well partition function and a model density of localised 
states, we find the temperature dependence of the specific 
heat. Under certain conditions the predicted specific heat 
follows ppecisely the sub-linear temperature dependence 
observedll3 in La2-xSrxCuOi at low temperatures. 



II. PARTITION FUNCTION 

The picture of interacting bosons with a short-range in- 
teraction filling up all the localised single-particle ptf-te 
in a random potential is known in the literatureEoilj. 
To calculate the partition function of localised bosons 
one has to take into account the repulsion between them. 
One cannot ignore the fact that the localisation length ^ 
generally varies with energy and diverges at the mobility 
edge. One would expect that the number of hard core 
bosons in a localised state near the mobility edge diverges 
in a similar way to the localisation length. Only a repul- 
sive interaction can stop all particles condensing into the 
lowest localised state. Thus, as stressed by Fisher et aZEj, 
for a Bose gas (in a random field) there is no sensible non- 
interacting starting point about which to perturb; this is 
in contrast to the case of a Fermi gas. 

Although the comprehensive scaling analysis of 
neutrafj and charged bosonsEil allows us to describe the 
quantum Bose glass-superfluid transition, the shape of 
the scaling functions as well as the thermodynamics of 
each phase away from the transition repiain unknown. 
A physically plausible model of neutralE£l and charged 
bosonaij in a random potential based on a separation 
of localised single-particle states frontjielocalised states 
might be helpful. As discussed in Ref.o, such separation, 
strictly speaking, is not possible as the localised states 
will be subsumed into the collective mode. Nevertheless, 



the effect of this mode can be understood by analysing 
the high density Umit where the dimensionless Coulomb 
repulsion is small: r^ = 4me^/eo(47rn/3)^/"^ < 1. Here m 
is the boson mass, n is the density and eo the static di- 
electric constant of the host material (and we take ?i = 1). 
In this limit thcpexcitation spectrum at T = was cal- 
culated by FoldytlLwho worked at zero temperature us- 
ing the Bogoliubovll3 approach. The Bogoliubov method 
leads to the result that the elementary excitations of the 
system have, for small momenta, energies characteristic 
of plasma oscillations which pass over smoothly for large 
momenta to the energies characteristic of single parti- 
cle excitations. They screen an external charge with a 
screening radius, JKhich is temperature independent at 
low temperaturestBEy. At high density, rs ^ 1, one can 
expect a disappearence of localised states because of the 
screening similar to the Mott metal-insulator transition 
in doped semiconductors. 

It is more likely, however, that the dimensionless 
strength of the Coulomb repulsion is about unity or 
larger, r<, > 1, and it is only when Vg is very large 
Tg >>> 1 that the Wigner crystallisation of charged 
carriers occurs. Hence, there is a wide interval of densi- 
ties where the localised states coexist with a superfluid. 
Being screened they depend on the interaction and the 
superfluid density. Nevertheless, based on the screening 
in the high density limittSEj we expect that the den- 
sity of localised states pL near the mobility edge, Ec 
remains temperature independent at low temperatures 
{ksTpi^ <C 1). If this is so, then the renormalisation 
of the effective single-particle energies by the collective 
mode does not affect the temperature dependence of any 
of the thermodynamic functions at low temperatures. 

With all these reservations we assume that at some fi- 
nite temperature Tc, bosons are condensed aX E — Ec so 
that the chemical potential ii = Ec- The excitation spec- 
trum of the delocalised superfluid has a gap of the order 
of the plasma frequencytll and so the Bogoliubov collec- 
tive modes can be ignored in the thermodynamics of the 
system (in iD their contribution is exponentially small 
while in 2D their energy scales as T^ and the specific heat 
as T"*). Even in the case of a short range repulsion the 
sound modes yield an energy proportional to T'^'^^ and 
hence a specific hea±_which behaves like C oc T'' (where 
d = dimensionality )Ej. Thus, for d > 2, the contribution 
to thermodynamics from the delocalised bosons appears 
to be negligible at low temperatures compared with that 
from bosons localised in shallow potential wells. Hence 
in the following we calculate the partition function and 
specific heat of localised bosons only. 

For simplicity, we choose 



C« 



Ec^Q 



(1) 



When two or more charged bosons are in a single localised 
state of energy E there may be significant Coulomb en- 
ergy and we try to take this into account as follows. The 
localisation length ^ is thought to depend on E via 



-Ef 



(2) 



where v > Q. The Coulomb potential energy of p charged 
bosons confined within a radius ^ can be expected to be 



potential energy 



P(p-l)e'^ 



(3) 



Thus the total energy of p bosons in a localised state of 
energy E is taken to be 



w{E) =pE + p{p - 1)k {-EY 



(4) 



where k > 0. Hence we see that the behaviour of charged 
bosons in localised states can be thought of as interme- 
diate between Bose-Einstein statistics and Fermi-Dirac 
statistics. When k = we have an equally spaced set 
of levels, i.e. Bose-Einstein behaviour, whereas when 
K = cxD we have Fermi-Dirac behaviour since the only 
levels with finite energy are p = and p — 1^ thus en- 
forcing an exclusion principle. When < k < oo we 
have the intermediate 'parastatistics' that the level spac- 
ing 5w — > oo as p — > oo. 

We take the total energy of charged bosons in localised 
states to be the sum of the energies of the bosons in 
the individual potential wells. The partition function Z 
for such a system is then the product of the partition 
functions for each of the wells, and the system free energy 
F = ksTluZ is simply the sum of the individual free 
energies kBT\Yiz{E). The free energy of the localised 
bosons in one unit cell is then given by 



F = kBT / dEpL (E) In z (E) , 



(5) 



where pL (E) is the one-particle density of localised states 
per unit cell below the mobility edge. 

In this paper we assume the following physically plau- 
sible density of localised states pl{E): 



Pl{E) 



7 



-e 1 



(6) 



We shall see below that, in the most interesting case, if 
the width of the tail 7 is large compared with fc^T, the 
specific heat is insensitive to the details of the shape of 
Pl{E), depending only on the value of pl{0)- 



III. THE THERMODYNAMICS OF A SINGLE 
POTENTIAL WELL 

Before considering the thermodynamics of charged 
bosons in our assumed density of localised states, we fir^ 
summarise here the results we have previously obtainedE£l 
for the thermodynamics of a single potential well with 
only one single particle level of energy e, i.e. for the case 



PLiE) = nLd{E - e) 



(7) 



The probability for this state to be occupied by p bosons 
is proportional to 

„-P{w{e)-p^i} 



where /x is the chemical potential and (3 
can re-express w — p/i as 



l/fcsT. We 
w-pn = K{-eY {p-paf - K{-ef pI (8) 



where 



1 



2K(-e)'- 



(9) 



Fig 1 shows a graph of w — pfi as a function of p. The 
partition function z{e) for such a single localised state is 



^w = E' 



-l3{w-pn) 
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pPo/3 
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(10) 



p=0 



The partition function is thus completely determined by 
the dimensionless parameters po a-nd ksT /[^{—eY]. 
The corresponding mean occupancy (jp) given by 



(p) = fcsT 



91nz(e) 
9/i 



and specific heat capacity c given by 







2 g^lnz(e) 
5/32 ■ 



(11) 



(12) 



are shown in Figs 2 and 3 for the case /i = 0. These 
results were calculated by simply truncating the partition 
function series at 100 terms. 

We now attempt to understand these results in more 
detail, looking separately at each temperature range. 

1. kBT<^n{-eY 

At low temperatures the partition function is dom- 
inated by the term with p closest to pq, i.e. the 
value of p giving the lowest value of w — pp., and 
so the mean occupancy {p) is an integer and goes 
up in steps as po increases, as seen in Fig 2. The 
changeover in dominance from one term to another 
occurs when po is a half-integer, at which point the 
two lowest energy states are degenerate. 

So long as one term dominates the partition func- 
tion, the specific heat c will be close to zero. How- 
ever when po is close to a half-integer we have a two 
level system and a corresponding Schottky anomaly 
in the specific heat capacity. Hence, at fixed tem- 
perature, the low temperature specific heat capac- 
ity (i.e. fc^T <C Pqk(— e)'^) is periodic in po. 



2. fcsT > K{-eY 

We can approximate the sum by an integral 



z{e) 



„p^/3K(-e)- 



dpe 



-/3K(-e)-(p-po)2 



(13) 



• n{-eY <kBT <plK{-eY 
In this case we can approximate the lower 
limit of the integral as — oo, i.e. the partition 
function can be approximated by an untrun- 
cated gaussian, and is therefore approximately 
symmetrical about po. Hence, in this temper- 
ature range we have 



{p) ~ Po 
as is clearly seen in Fig 2, and 



C^\kB 



(14) 



(15) 



as seen in Fig 3. 

kuT > K{-eY and kgT » p^ni-eY 
In this case we can approximate the partition 
function as being an integral over half a Gaus- 
sian. We then obtain 



ip) 



I ksT 



(16) 



and the specific heat is again given by c w 
ks/^; both of these results are seen in Figs 2 
and 3. 



IV. THE NUMBER OF BOSONS IN LOCALISED 

STATES 

Having established the partition function z{E) for a 
single potential well containing one single particle level 
of energy E we now use equations and o to derive the 
thermodynamics of localised bosons. 

The average number Nl of bosons in localised states 
in each unit cell is 



Nr 



dE{p)pL{E) 



(p) can be approximated as 




if keT < pIk{-EY 



— 7^-w^ otherwise 



(17) 



(18) 



Hence a reasonable approximation for (p) at all temper- 
atures is 



{p)^l + 



1 



knT 



2 2k{~EY-^ y t:k{~EY 
With this approximation we obtain 



(19) 



(20) 



Hence if z^ > 2 then iV^ becomes infinite and so the 
formation of a superfluid is excluded. 



V. THE SPECIFIC HEAT OF LOCALISED 
BOSONS 

We have shown above that the contribution to the spe- 
cific heat c from a single localised state is determined 
solely by the values of two dimensionless quantities: Po 
and ksT / [k,{—EY\. If v is temperature independent, 
c can instead be considered as a function of a different 
pair of quantities: —E/ [ksT] and t = k {ksTY^ The 
integrated specific heat C of the localised bosons is then 



C 



dEpL{EY 



-E 
k^' 



(21) 



A. Specific heat for v = 1 

In this case t = k, and is therefore temperature inde- 
pendent. Thus we obtain: 



C^ksT 



-E 



knT 



pl{ey 



-E 



knT' 



(22) 



If pl{E) is constant in the region —E <r-^ ksT/n then 
equation ^ yields C (x T; for p^ given by equation |g this 
result applies if 7 » ksT / k. 



This region only exists if r < 1. c is dominated by 
Schottky anomalies, each centred on a half-integer 
value of Po. The last full Schottky anomaly in the 
series is at po = | , so the total number s of Schot- 
tky anomalies is 



Vi/^_i 



(24) 



If we ignore all except the lowest two energy levels 
then for each anomaly 



ibove 



dEc 



below 



-E 
k^' 



3(^-1) 



klT 



(25) 



3. -E/ [keT] > T 



-1/(^-1) 



{i.e. Po < 1 } 



Here the p = 1 level is the lowest energy level, and 
the p = level is the next lowest. If higher levels are 
neglected then c is solely a function of —E/ [fc^T] 
and is non-negligible only for 0.1 < —E/ [fc^T] < 
10. Hence this region can be neglected entirely if 

T < (jq)"' ■ If, on the other hand, r > 10""^ the 
integrated specific heat for this region will be that 
of half a Schottky anomaly: 



J regio7i 3 \ '^B ^ 



1 



■kiT 



(26) 



If K'y'^~^ ^ 1 then e^^"* ~ 1 throughout regions (1) 
and (2). If 7 :^ ksT then e^/'*' ~ 1 throughout region 
(3). We take both of these conditions to hold so that 
equation El^ can be re- written as 



C 
T 



riLki 

1 



-E 



knT 



-E 



(27) 



i.e. if V is fixed then C /T is solely a function of r. 

Hence we can distinguish two cases for which an ap- 
proximate analytical formula for C can be derived: 



B. Specific heat for v > 1 

For v > 1 the integral in equation |l] can be thought 
of as an integral over three regions: 

1. -E/ [keT] < r-i/- {i.e. p„ > i (l -^ r-V-) } 
In this region fcsT > k{—EY and so 



— kn 
2 ^ 



(23) 



1. T< 






In this case region (3) can be neglected. From equa- 
tions P3, H and E3 we then obtain 



(28) 
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r ^2 ^1 


nLk% 
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tI/i'7 


This implies 







C oc T^/"" 



(29) 



2. T-i/'' < -E/ [keT] < r-i/(— 1) 

{i.e. ^{l + T-^/") >Po>l} 



i.e. the specific heat has a power law dependence on 
temperature in which the power is less than unity. 



2. T> 10"-^ 

In this case the specific heat C is dominated by 
region (3). From equation Era we then obtain 



T 



67 



(30) 



i.e. C /T is temperature independent. 



Fig. 4 compares these two formulae (for the case v = 
1.5) with the result of calculating C from equation E7^ by 
numerical integration (in which c is calculated from the 
first 100 terms of the partition function). 



C. Specific heat for 1/ <1 

Now Po increases as — -E increases. Again we can dis- 
tinguish three regions 

1. -E/ [ksT] < r-i/- {i.e. p^ < 1 (1 + r-i/^) } 
Here ksT > k{—EY and so once again c « kB/2. 

2. T-i/'' < -£;/ [fcsT] < ri/(i"'') 

{i.e. i(l + T-i/'') <Po<l} 

Here the partition function is dominated by the two 
levels corresponding to p = 1 and p = 0. If t < 1 
this region does not exist. 

3. -E/ [keT] > ri/fi-") {i.e. Po < I } 

Here c is dominated by an infinite series of Schottky 
anomalies centred on values of E corresponding to 
half-integer values of po^ 



E = [i2po - 1)k,] 



(31) 



We consider the low temperature limit, i.e. t ^ 1 (for 
v < 1). Let 



X 



-E 

7 



Equation G^ becomes 



C = m 



'^cdx 



(32) 



(33) 



In the low temperature limit ksT ^ i^l'^ , region 1 makes 
a negligible contribution to the integral in equation p3. 
Provided ksT <C (1 — 1^)7 we can make the approxi- 
mation that, when integrating through each Schottky 
anomaly, e~^ can be taken to be a constant depending on 
the anomaly; hence the contribution from each anomaly 
is proportional to the result given in equation Gq (and 
equation E9). Thus, in the low temperature limit, we 
have C (xT. 



If 7^""^ ^ K then many Schottky anomalies contribute 
significantly to the integral in equation ^. In that case 
we can neglect region 2 and transform the sum over 
Schottky anomalies into an integral with respect to po, 
obtaining: 



C 



3(^ - 1)7 Jo 

6^7*^ 
Tr^klnLTjl 
6K7'' 



^dpo 



X ''e ^dx 



-T 



(34) 



In the high temperature limit, on the other hand, we 
have T ^ 1 and ksT ':^ k^'^; in that case only region 1 
is significant and we obtain: 

C « -ULkB (35) 



VI. CONCLUSION 

We have used a reasonable scaling of Coulomb energy 
with localisation length to calculate the low temperature 
dependence of the specific heat of a charged superfluid in 
a random potential. The result strongly depends on the 
exponent i' of the localisation length. While the specific 
heat is linear for 1/ < 1, it is proportional to T^^'^ for 
i^ > 1. 

We believe that our findings are relevent for those 
doped high- Tc cupjcates having many properties of a 
charged Bose-liquidtl. Because of the high level of doping 
of these Mott-Hubbard insulators in the superconduct- 
ing region one can expect a continuous density of states 
Pl{E) similar to that studied by us; our prediction about 
the power law dependence of C (for v > \) requires, of 
the density of states, only that pl{Ec) ^0. j— . 

Superconducting La2-xSrxCuOi has been observedll2l 
to have a low temperature specific heat proportional to 
T" with a being in the range 0.33 < a < 0.78 depen- 
dent on doping. This fits with our prediction for charged 
bosons partly localised by disorder at low temperature 
and V > 1. The doping dependence of a — Ijv could 
be due to a gradual change of the localisation length ex- 
ponent V caused by a difference in the random potential 
profiles for different x. |— , 

Values of a as low as 0.33 are observedty experimen- 
tally while our model excludes the possibility of super- 
conductivity for a = X/v < 0.5. However we believe 
that the approximation for Coulomb potential energy in 
equations 0-0 is likely to be inadequate at large ^. 



Enlightening discussions with N. Hussey, V. Kabanov, 
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FIG. 1. Graph of {w ~ p^)/[n{-EY] against p. 
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FIG. 2. The mean occupancy (p) of a single level as a func- 
tion of po and logj^g {fcflr/[K(— e)"]}. 



FIG. 3. The specific heat capacity c of a single level as a 
function of po and logjQ {fc_Br/[fi:(— e)"]}. 



FIG. 4. Specific heat capacity C divided by temperature T 
for the case f = 1.5. The continuous curve is obtained from 
numerical integration of the integral in equation ETI c being 
evaluated by direct summation of the first 100 terms of the 
partition function. The two dashed lines are the predictions 
from equations Eq and pCj. 
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